Systematic alteration of ATAC-seq for profiling open chromatin in cryopreserved nuclei preparations from livestock tissues.

The use of Assay for Transposase-Accessible Chromatin (ATAC-seq) to profile chromatin accessibility has surged over the past years, but its applicability to tissues has been very limited. With the intent of preserving nuclear architecture during long-term storage, cryopreserved nuclei preparations from chicken lung were used to optimize ATAC-seq. Sequencing data were compared with existing DNase-seq, ChIP-seq, and RNA-seq data to evaluate library quality, ultimately resulting in a modified ATAC-seq method capable of generating high quality chromatin accessibility data from cryopreserved nuclei preparations. Using this method, nucleosome-free regions (NFR) identified in chicken lung overlapped half of DNase-I hypersensitive sites, coincided with active histone modifications, and specifically marked actively expressed genes. Notably, sequencing only the subnucleosomal fraction dramatically improved signal, while separation of subnucleosomal reads post-sequencing did not improve signal or peak calling. The broader applicability of this modified ATAC-seq technique was tested using cryopreserved nuclei preparations from pig tissues, resulting in NFR that were highly consistent among biological replicates. Furthermore, tissue-specific NFR were enriched for binding motifs of transcription factors related to tissue-specific functions, and marked genes functionally enriched for tissue-specific processes. Overall, these results provide insights into the optimization of ATAC-seq and a platform for profiling open chromatin in animal tissues.

Mapping of the chicken cleft primary palate mutation on chromosome 11 and sequencing of the 4.9 Mb linked region.

An embryonic lethal mutation in chicken named cleft primary palate (cpp) is inherited in an autosomal recessive mode and results in a severely truncated upper beak. In this study, genotyping and sequencing techniques were employed to advance our genetic and genomic knowledge of the mutation's chromosomal location, candidate region and possible causative element using a congenic inbred line. Herein, the candidate region for the cpp developmental mutation was established as a ca. 5.1 Mb region of chicken chromosome 11 (GGA 11) through the use of a 600K Affymetrix SNP array. The SNPs identified from this array linked to cpp were used to genotype individuals from the congenic inbred line over several generations and thereby fine-map the causative region resulting in an approximately 200 kb size reduction. This candidate region (4.9 Mb) was sequenced via capture array in a cohort of 24 individuals, including carriers, mutants and their wild type (wt) siblings. Interestingly, the GGA 11 region for cpp encompasses the predicted centromere location and is thus unlikely to be highly disrupted by further recombination. Here we report on the variation unique to the cpp mutation, i.e. single-nucleotide variants and insertions or deletions. Although the candidate region contains several genes of interest with regard to the cpp phenotype, only one cpp-linked variant was predicted to have a significant physiological effect by causing a frameshift mutation in ESRP2, which has a role in tissue-specific splicing during development.

Narrowing the wingless-2 mutation to a 227 kb candidate region on chicken chromosome 12.

Wingless-2 (wg-2) is an autosomal recessive mutation in chicken that results in an embryonic lethal condition. Affected individuals exhibit a multisystem syndrome characterized by absent wings, truncated legs, and craniofacial, kidney, and feather malformations. Previously, work focused on phenotype description, establishing the autosomal recessive pattern of Mendelian inheritance and placing the mutation on an inbred genetic background to create the congenic line UCD Wingless-2.331. The research described in this paper employed the complementary tools of breeding, genetics, and genomics to map the chromosomal location of the mutation and successively narrow the size of the region for analysis of the causative element. Specifically, the wg-2 mutation was initially mapped to a 7 Mb region of chromosome 12 using an Illumina 3 K SNP array. Subsequent SNP genotyping and exon sequencing combined with analysis from improved genome assemblies narrowed the region of interest to a maximum size of 227 kb. Within this region, 3 validated and 3 predicted candidate genes are found, and these are described. The wg-2 mutation is a valuable resource to contribute to an improved understanding of the developmental pathways involved in chicken and avian limb development as well as serving as a model for human development, as the resulting syndrome shares features with human congenital disorders.

Virus and host genomic, molecular, and cellular interactions during Marek's disease pathogenesis and oncogenesis.

Marek's Disease Virus (MDV) is a chicken alphaherpesvirus that causes paralysis, chronic wasting, blindness, and fatal lymphoma development in infected, susceptible host birds. This disease and its protective vaccines are highly relevant research targets, given their enormous impact within the poultry industry. Further, Marek's disease (MD) serves as a valuable model for the investigation of oncogenic viruses and herpesvirus patterns of viral latency and persistence--as pertinent to human health as to poultry health. The objectives of this article are to review MDV interactions with its host from a variety of genomic, molecular, and cellular perspectives. In particular, we focus on cytogenetic studies, which precisely assess the physical status of the MDV genome in the context of the chicken host genome. Combined, the cytogenetic and genomic research indicates that MDV-host genome interactions, specifically integration of the virus into the host telomeres, is a key feature of the virus life cycle, contributing to the viral achievement of latency, transformation, and reactivation of lytic replication. We present a model that outlines the variety of virus-host interactions, at the multiple levels, and with regard to the disease states.

The expression of preaxial polydactyly is influenced by modifying genetic elements and is not maintained by chromosomal inversion in an avian biomedical model.

Polydactyly (Po) is a common mutation found in many vertebrates. The UCD-Po.003 congenic chicken line was previously characterized for Po inheritance (autosomal dominant) and the mutation was mapped to chromosome 2p. Here, we describe for the first time the range and variability of the phenotype in this congenic line. Further, we studied the hypothesis that a chromosomal inversion was responsible for the maintenance of a large (6.3 Mb) candidate gene region. Fluorescence in situ hybridization employing BACs encompassing a 10.7-Mb region of GGA2p showed that the Po chromosome was normal, i.e. exhibits the wild-type BAC order. Continued fine-mapping along with a change in breeding strategy reduced the size of the causative region to 1.43 Mb. Recent research indicates that the cause of preaxial Po resides within a 794-bp highly conserved zone of polarizing activity regulatory sequence (ZRS) element located in intron 5 of the LMBR1 gene; however, the ZRS polymorphism of interest is found in some but not all breeds of polydactylous chicken. Therefore, we sequenced the ZRS in 101 heterozygous and 30 unaffected (wild-type) individuals to establish the relevance of this region to the Po condition in the UCD-Po.003 congenic line. A single point mutation (C/A at coordinate GGA2p: 8,414,121) within the ZRS segregated with carrier status. The polydactylous UCD-Silkie line also maintains this SNP in addition to a single base deletion. An inheritance analysis of the phenotypic variation in UCD-Po.003 suggests recessive epistasis as the mode of inheritance for the additional modifying genetic elements, residing outside the ZRS, to impact the preaxial polydactyl phenotype. These results contribute to our understanding of the cause of Po in an important vertebrate model.

Molecular and cellular evidence for the alternative lengthening of telomeres (ALT) mechanism in chicken.

Telomere maintenance is an important genetic mechanism controlling cellular proliferation. Normally, telomeres are maintained by telomerase which is downregulated upon cellular differentiation in most somatic cell lineages. Telomerase activity is upregulated in immortalized cells and cancers to support an infinite lifespan and uncontrolled cell growth; however, some immortalized and transformed cells lack telomerase activity. Telomerase-negative tumors and immortalized cells utilize an alternative mechanism for maintaining telomeres termed alternative lengthening of telomeres (ALT). This research explored evidence for the ALT pathway in chicken cell lines by studying nontransformed immortalized cell lines (DF-1 and OU2) and comparing them to a normal (mortal) cell line and a transformed cell line (DT40). The research consisted of molecular and cellular analyses including profiling of telomeric DNA (array sizing and total content), telomerase activity, and expression of genes involved in the telomerase, recombination, and ALT pathways. In addition, an immunofluorescence analysis for an ALT marker, i.e. ALT-associated promyelocytic leukemia bodies (APBs), was conducted. Evidence for ALT was observed in the telomerase-negative immortalized cell lines. Additionally, the APB marker was also found in the other cell systems. The attributes of the chicken provide an additional vertebrate model for investigation of the ALT pathway.

Poultry genome sequences: progress and outstanding challenges.

The first build of the chicken genome sequence appeared in March, 2004 - the first genome sequence of any animal agriculture species. That sequence was done primarily by whole genome shotgun Sanger sequencing, along with the use of an extensive BAC contig-based physical map to assemble the sequence contigs and scaffolds and align them to the known chicken chromosomes and linkage groups. Subsequent sequencing and mapping efforts have improved upon that first build, and efforts continue in search of missing and/or unassembled sequence, primarily on the smaller microchromosomes and the sex chromosomes. In the past year, a draft turkey genome sequence of similar quality has been obtained at a much lower cost primarily due to the development of 'next-generation' sequencing techniques. However, assembly and alignment of the sequence contigs and scaffolds still depended on a detailed BAC contig map of the turkey genome that also utilized comparison to the existing chicken sequence. These 2 land fowl (Galliformes) genomes show a remarkable level of similarity, despite an estimated 30-40 million years of separate evolution since their last common ancestor. Among the advantages offered by these sequences are routine re-sequencing of commercial and research lines to identify the genetic correlates of phenotypic change (for example, selective sweeps), a much improved understanding of poultry diversity and linkage disequilibrium, and access to high-density SNP typing and association analysis, detailed transcriptomic and proteomic studies, and the use of genome-wide marker- assisted selection to enhance genetic gain in commercial stocks.

Chromosomal mapping of chicken mega-telomere arrays to GGA9, 16, 28 and W using a cytogenomic approach.

Four mega-telomere loci were mapped to chicken chromosomes 9, 16, 28, and the W sex chromosome by dual-color fluorescence in situ hybridization using a telomeric sequence probe and BAC clones previously assigned to chicken chromosomes. The in-common features of the mega-telomere chromosomes are that microchromosomes are involved rather than macrochromosomes; in three cases (9, 16, 28) acrocentrics are involved with the mega-telomeres mapping to the p arms. Three of the four chromosomes (9, 16, W) encode tandem repeats which in two cases (9 and 16) involve the ribosomal DNA arrays (the 5S and 18S-5.8S-28S gene repeats, respectively). All involved chromosomes have a typical-sized telomere on the opposite terminus. Intra- and interindividual variation for mega-telomere distribution are discussed in terms of karyotype abnormalities and the potential for mitotic instability of some telomeres. The diversity and distribution of telomere array quantity in the chicken genome should be useful in contributing to research related to telomere length regulation - how and by what mechanism genomes and individual chromosomes establish and maintain distinct sets of telomere array sizes, as well as for future studies related to stability of the chicken genome affecting development, growth, cellular lifespan and disease. An additional impact of this study includes the listing of BAC clones (26 autosomal and six W BACs tested) that were cytogenetically verified; this set of BACs provide a useful tool for future cytogenetic analyses of the microchromosomes.

Avian genetic stocks: the high and low points from an academia researcher.

Poultry genetic resources, which are valuable for research, span an impressive gamut from breeds to highly specialized inbred lines. The community of scientists utilizing specialized lines is broad, including researchers in medicine, basic biology, and agricultural science. The majority of specialized research lines used by such scientists are held at land grant universities. Over the prior 2 decades, hundreds of lines were eliminated. This pattern continues today with no evidence of abatement. Awareness and visibility of the causes and ongoing problems have been highlighted via a number of high-profile forums. Given the large community of scientists and the negative impact on future advances in biological, medical, and agricultural research as these genetic resources dwindle, the issue is of national interest and warrants federal funding to support a united network of avian and poultry stocks centers.

Meiotic instability of chicken ultra-long telomeres and mapping of a 2.8 megabase array to the W-sex chromosome.

The objective of this research was to study the meiotic stability of a subset of chicken telomere arrays, which are the largest reported for any vertebrate species. Inheritance of these ultra-long telomere arrays (200 kb to 3 mb) was studied in a highly homozygous inbred line, UCD 003 (F >or= 99.9). Analysis of array transmission in four families indicated unexpected heterogeneity and non-Mendelian segregation including high-frequency-generation of novel arrays. Additionally, the largest array detected (2.8 Mb) was female-specific and correlated to the most intense telomeric DNA signal on the W-sex chromosome by fluorescence in situ hybridization (FISH). These results are discussed in regard to the potential functions of the ultra-long telomere arrays in the chicken genome including generation of genetic variation through enhanced recombination, protection against erosion by providing a buffer for gene-dense regions, and sex-chromosome organization.

Telomeres in the chicken: genome stability and chromosome ends.

Telomeres are the complex nucleoprotein structures at the termini of linear chromosomes. Telomeric DNA consists of a highly conserved hexanucleotide arranged in tandem repeats. Telomerase, a ribonucleoprotein of the reverse transcriptase family, specifies the sequence of telomeric DNA and maintains telomere array length. Numerous studies in model organisms established the significance of telomere structure and function in regulating genome stability, cellular aging, and oncogenesis. Our overall research objectives are to understand the organization of the telomere arrays in chicken in the context of the unusual organization and specialized features of this higher vertebrate genome (which include a compact genome, numerous microchromosomes, and high recombination rate) and to elucidate the role telomeres play in genome stability impacting cell function and life span. Recent studies found that the chicken genome contains three overlapping size classes of telomere arrays that differ in location and age-related stability: Class I 0.5 to 10 kb, Class II 10 to 40 kb, and Class III 40 kb to 2 Mb. Some notable features of chicken telomere biology are that the chicken genome contains ten times more telomeric DNA than the human genome and the Class III telomere arrays are the largest described for any vertebrate species. In vivo, chicken telomeres (Class II) shorten in an age-related fashion and telomerase activity is high in early stage embryos and developing organs but down-regulates during late embryogenesis or postnatally in most somatic tissues. In vitro, chicken cells down-regulate telomerase activity unless transformed. Knowledge of chicken telomere biology contributes information relevant to present and future biotechnology applications of chickens in vivo and chicken cells in vitro.

The chicken telomerase RNA gene: conservation of sequence, regulatory elements and synteny among viral, avian and mammalian genomes.

Telomerase RNA (TR) is essential for telomerase activity and the maintenance of telomere length in proliferating cell populations. The objective of the present research was to define the cytogenetic and molecular genomic organization of chicken TR (chTR). The chTR exists as a single copy gene (TERC, alias TR), mapping to chromosome 9 (GGA9). The loci on the q arm of GGA9 map to three chromosomes in human with five of the nine GGA9q loci mapping to HSA3q. Sequencing of the chTERC locus (3,763 bp) from the UCD 001 genome (Red Jungle Fowl) included: 604 bp 5', 465 coding, and 2,694 bp 3' (from -604 to +3159). Sequence analysis included homology searches conducted on several levels including comparisons among different chicken genotypes, Marek's disease virus (MDV) sequences, plus human and murine. We provide evidence for distal 5' and 3' sequence homology between chTERC and the MDV genome among other known regions of homology (promoter and coding), elaborate on 5' transcription factor binding motifs among the various genomes as well as show type and number of TERT-related motifs 3' of chicken TR (e.g., Sp1, c-Myb, c-Myc, AP2, among others). Surrounding the gene are more than 25 Sp1 sites, over 20 oncogene transcription factor binding motifs and numerous hormonal and other specialized binding motifs. Knowledge of 5' and 3' chTERC regulatory elements will be useful for investigating normal control mechanisms during growth and development as well as investigating the potential for dysregulation of this important gene during oncogenesis, especially among different genotypes.

Dynamics of telomere erosion in transformed and non-transformed avian cells in vitro.

Although vertebrate telomeres are highly conserved, telomere dynamics and telomerase profiles vary among species. The objective of the present study was to examine telomerase activity and telomere length profiles of transformed and non-transformed avian cells in vitro. Non-transformed chicken embryo fibroblasts (CEFs) showed little or no telomerase activity from the earliest passages through senescence. Unexpectedly, a single culture of particularly long-lived senescent CEFs showed telomerase activity after over 250 days in culture. Transformed avian lines (six chicken, two quail and one turkey) and tumor samples (two chicken) exhibited telomerase activity. Telomere length profiles of non-transformed CEF cultures derived from individual embryos of an inbred line (UCD 003) exhibited cycles of shortening and lengthening with a substantial net loss of telomeric DNA by senescence. The telomere length profiles of several transformed cell lines resembled telomere length profiles of senescent CEFs in that they exhibited little of the typical smear of terminal restriction fragments (TRFs) suggesting that these transformed cells may possess a reduced amount of telomeric DNA. These results show that avian telomerase activity profiles are consistent with the telomerase activity profiles of human primary and transformed cells. Further, monitoring of telomere lengths of primary cells provides evidence for a dynamic series of changes over the lifespan of any specific cell culture ultimately resulting in net telomeric DNA loss by senescence.

Organization of telomere sequences in birds: evidence for arrays of extreme length and for in vivo shortening.

Telomeres are the specialized ends of chromosomes consisting of highly conserved repeat (5'-TTAGGG-3')(n) sequences. Lack of information regarding the existence of an in vivo telomere clock function in birds, conflicting data regarding telomere array length in the chicken model, and the paucity of molecular telomere information for other avian species led us to study telomere array organization within and among 18 species and subspecies of birds. Most of the species contained between 2% and 4% telomere sequence per diploid genome. Arrays spanning 0.5-10 kb (Class I) and 10-40 kb (Class II) were observed in all of the species studied. Extremely long arrays, ranging from hundreds of kilobases to 1-2 Mb (Class III) were observed in all except two raptor species, the northern goshawk and American bald eagle. In chicken, there was evidence for shortening of the Class II arrays in vivo, based on intraindividual comparisons of somatic versus germline tissues in birds of different ages; terminally differentiated erythrocyte arrays were, on average, 2.3 kb shorter than sperm (germline) arrays. This study provides the first evidence for the existence of telomere arrays significantly larger than have been described for any vertebrate species to date and for developmentally programmed in vivo telomere shortening in the Aves taxa. The novel finding of megabase-sized telomere arrays may be an important feature of avian karyotypes that contain a large number of very small genetic units, the microchromosomes.

Regulation of brain-derived neurotrophic factor messenger RNA levels in avian hypothalamic slice cultures.

Mechanisms regulating the expression of brain-derived neurotrophic factor, a member of the neurotrophin family, have been extensively studied in the rat cerebral cortex, hippocampus and cerebellum. In contrast, little is known regarding the regulation of this growth factor in the hypothalamus. Here we present an analysis of the regulation of brain-derived neurotrophic factor messenger RNA levels in chick embryo hypothalamic slice cultures following exposure to potassium chloride, glutamate agonists and sex steroids. Following a week in chemically-defined media the tissue was depolarized by exposure to 50 mM potassium chloride for 6h, resulting in a significant 4.2-fold increase in the level of brain-derived neurotrophic factor messenger RNA. This result is consistent with studies of other brain regions. Similar 6-h acute exposures of the hypothalamic cultures to 25 microM N-methyl-D-aspartic acid, 25 microM kainic acid and 25 microM alpha-amino-3-hydroxy-5-methyl-isoxazole-4-propionic acid also significantly increased messenger RNA levels 2.5-, 2.1- and 1.4-fold, respectively. It was previously reported that brain-derived neurotrophic factor levels within the rat cerebral cortex, olfactory bulb and hippocampus are altered by exposure to 17beta-estradiol. Here we show that in hypothalamic slice cultures neither acute nor chronic treatments with 10 and 100 nM 17beta-estradiol and 10nM testosterone significantly altered the steady-state level of this growth factor. These findings show that neuronal activity, induced by glutamate agonists and potassium chloride, can regulate brain-derived neurotrophic factor messenger RNA levels within embryonic hypothalamic slice cultures. This regulation could play a critical role in the modulation of programmed cell death and synaptic maturation during development of the hypothalamus.

Patterns of ribosomal gene variation in elite commercial chicken pure line populations.

The nucleolus organizer region (NOR) encodes the tandemly repeated 18S, 5.8S and 28S ribosomal (r) RNA genes. The NORs of broiler and layer commercial chicken pure lines were studied to establish the type and extent of genetic variation at this important locus. The parameters studied were gene copy number, repeat size, and diversity of NOR-types. The populations were organized into three groups for analysis including brown-egg broiler (13 lines), brown-egg layer (six lines), and white-egg layer (eight lines). The ribosomal gene copy number average of the white-egg layer populations was significantly lower (329 genes) than that of the brown-egg layers (372 genes); the brown-egg broiler ribosomal gene average was intermediate (350 genes). The white-egg layer populations exhibited a ribosomal repeat unit average size of 36 kb, significantly different from the brown-egg layer and brown-egg broiler average repeat unit size of 32.5 and 33.9 kb, respectively. NOR array size was similar among the three groups (6 mb). The brown-egg broiler populations exhibited polymorphic NOR patterns, intra- and interline, whereas the white-egg layer populations were essentially monomorphic for NOR-type; brown-egg layers exhibited an intermediate level of NOR diversity. Some NOR array characteristics may be a function of breed origin as brown-egg commercial populations, both broilers and layers, have similar breed origins and exhibited similarities for predominant repeat unit size as compared with white-egg layer populations. However, the finding that brown-egg broiler lines typically exhibit a greater number of segregating NOR-types than brown-egg layer lines suggests that the selection schemes of broiler vs. layer pure line populations may also have influenced the degree of variation at this gene complex.

Ontogeny of telomerase in chicken: impact of downregulation on pre- and postnatal telomere length in vivo.

Telomeres are the termini of linear chromosomes composed of tandem repeats of a conserved DNA sequence. Telomerase provides a mechanism for proliferating cells to offset telomeric sequence erosion by synthesizing new repeats onto the end of each parental DNA strand. Reduced or absent telomerase activity can lead to telomere shortening and genome instability. Telomeres and telomerase have not previously been characterized during ontogeny of any avian species. In the present study, telomerase activity in the chicken model was examined from early differentiation embryos through to adulthood. Telomerase activity was detected in all early embryos (preblastula through neurula) and in tissues throughout organogenesis. Subsequently, telomerase was downregulated in the majority of somatic tissues, either pre- or postnatally. A subset of tissues, such as intestine, immune and reproductive organs, exhibited constitutive activity. The impact of telomerase downregulation on telomere length was investigated and a telomere reduction of 3.2 kb in somatic tissues compared with germ line was observed in 5-year-old adults. The present results suggest that the telomere clock function is a conserved feature of avians as well as mammals. Knowledge regarding the relationships among telomerase regulation, proliferation/senescence profiles and differentiation status will be useful for numerous applications of chicken cells.

Molecular characterization of ribosomal gene variation within and among NORs segregating in specialized populations of chicken.

The molecular organization of the 18S, 5.8S, and 28S ribosomal RNA gene repeat units, located at the single nucleolus organizer region (NOR) locus in the chicken, was investigated in genetically distinct populations of research and commercial chickens. Substantial gene repeat variation within and among NORs was documented. Intact ribosomal gene repeat size ranged from 11 kb to over 50 kb. Unique combinations of ribosomal genes, of different size, were specific to particular populations. It was determined that the basis for the ribosomal gene repeat size variation was intergenic spacer (IGS) length heterogeneity. Interestingly, in different populations, the location of the variation that contributes to length heterogeneity was specific to particular IGS subregions. In addition to IGS variation, an inbred line of Red Jungle Fowl exhibited coding region variation. Ribosomal gene copy number variation was also studied, and line averages ranged from 279 to 368. Average rDNA array size (a function of copy number and gene repeat length) was calculated for each of the populations and found to vary over a range of two megabases, from 5 to 7 Mb.